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Abstract ..."but we do not have quantum gravity." This phrase is often used when analysis of a 
physical problem enters the regime in which quantum gravity effects should be taken into account. 
In fact, there are several models of the gravitational field coupled to (scalar) fields for which the 
quantization procedure can be completed using loop quantum gravity techniques. The model we 
present in this paper consist of the gravitational field coupled to a scalar field. The result has similar 
| structure to the loop quantum cosmology models, except for that it involves all the local degrees of 

freedom of the gravitational field because no symmetry reduction has been performed at the classical 
level. 
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O ! I. INTRODUCTION 

CT 

H . The recent advances in loop quantum gravity (LQG) UQ, strongly suggest that the goal of 

constructing a candidate for quantum theory of gravity and the Standard Model is within reach. 
Remarkably, that goal can be addressed within the canonical formulation of the original Einstein's 
general relativity in four dimensional spacetime. A way to define 'physical' dynamics in a back- 
ground independent theory, where spacetime diffeomorphisms are treated as a gauge symmetry, is 
the framework of relational Dirac observables (often also called "partial" observables 0,0, QiB 
section 1.2 of Q). The main idea is, that part of the fields adopt the role of a dynamically coupled 
. observer, with respect to which the physics of the remaining degrees of freedom in the system is 

formulated. In this framework the emergence of the dynamics, time and space can be explained as 
an effect of the relations between the fields. As far as technical issues of a corresponding quantum 
' theory are concerned, the most powerful example of the relational observables framework is the 

deparametrization technique [9rll2j . This allows to map canonical General Relativity into a theory 
with a (true) non-vanishing Hamiltonian, that is independent of the (emergent) time provided by 
the observer fields. All this can be achieved at the classical level, the framework of Loop Quan- 
tum Gravity (LQG) itself, provides then the tools of the quantum theory like quantum states, 
the Hilbcrt spaces, quantum operators of the geometry and fields and well defined quantum op- 
erators for the classical constraints of General Relativity (see 0,(3] and references therein). The 
combination of LQG with the relational observables and deparametrization framework makes it 
possible to construct general relativistic quantum models. Applying LQG techniques to perform 
the quantization step has the consequence that the quantum fields of the Standard Model have 
to be reintroduced within the scheme of LQG. This is due to the reason that the standard quan- 
tum field theory (QFT) defined on the Minkowski (or even ADS) background is incompatible with 
quantization approach used in LQG. Therefore, the resulting quantum theory of gravity cannot 
be just coupled to the Standard Model in it's present form. The formulation of the full Standard 
Model within LQG will require some work. For this reason, we proceed step by step, increasing 
gradually the level of complexity. The first step was constructing various cosmological models by 
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analogy with LQG by performi ng a sy mmetry reduction already at the classical level. They give 
rise to loop quantum cosmology [LML9I ] (LQC). We have learned from them a lot about qualitative 
properties of quantum spacetime and its quantum dynamics [20I . [2~l| . That knowledge is very useful 
in performing the second step, that is introducing quantum models with the full set of the local 
gravitational degrees of freedom. The first quantum model of the full, four dimensional theory 
of gravity was obtained by applying LQG techniques [23] to the Brown-Kuchar model of gravity 
coupled to dust Q. In the current paper we apply LQG to the model introduced by Rovelli and 
Smolin [24| whose classical canonical structure was studied in detail by Kuchar-Romano [25[ . This 
is a model of gravity coupled to a massless scalar field. Our goal is to complete the construction of 
the quantum model with the tools of LQG. In the firs part of the paper (Sections I-II) we introduce 
the model, study the structure of the space of solutions to the quantum constraints, and the Dirac 
observables, assuming only suitable Hilbert products and operators exist. The result of this part is 
a list of mathematical elements necessary and sufficient for the model to exist. In the second part 
(Section III) we apply the framework of LQG. We show it provides the necessary Hilbert spaces 
and operators, and complete the construction of the model. 



II. CANONICAL GRAVITY COUPLED TO A CLASSICAL SCALAR FIELD 



A. The standard approach 



The point of our interest in this paper is gravity coupled to a scalar field. We are considering a 
metric tensor field q a b and a scalar field cf> on a 3-manifold M (the space). The conjugate momenta 
are denoted respectively by p ab and ir. The only non- vanishing Poisson brackets among them are 

{q ab (x),p cd iy)} = 5{x,y)8lX), {<f>(x)Mv)} = (2-1) 

The intrinsic and extrinsic geometry of M (as M being the Cauchy surface of 4-dimcnsional space- 
time) is described by the first pair of canonically conjugate variables (q a b,p ab )- The field q a b defines 
the intrinsic Riemann geometry of M whereas p ab contains the information about the extrinsic cur- 
vature of M imbedded in the spacetime. 

The variables (q a b,p ab ) are known from the standard canonical formulation of gravity usually 



called ADM formalism 33[ (see also chapter 10 and appendix E of [26|). But one can use any other 
variables in this part of our paper (Section I-II). In Section III we will apply loop quantum gravity 
(LQG), and therein we will be using the Ashtckar-Barbcro variables (A % a ,P"), i = 1, 2, 3 (and the 
notation of Q). They are also canonically conjugate to each other, and the only non- vanishing 
Poisson bracket is 

{Ai(x), P b (y)} = 5(x,y)6 b a 6>. (2.2) 

The intrinsic and extrinsic geometry of M can be recovered out of them, as they are defined by the 
orthonormal coframe e^, the corresponding connection 1-form T l . the extrinsic curvature 1-form 
K l a and a fixed Barbero-Immirzi parameter 7 (for its value see |27H29j | ) . namely 

Al = r a+1 K l a , P? = —^-eieU bc ^k (2.3) 

where 7/ 123 = 1 = 6123 and r] abc , e a bc are completely antisymmetric. 

The fields (A l a ,P°-) set an su(2) valued 1-form, and, respectively, su(2)* valued vector density 

A= A\{x)n®dx a , P = P?(x)t 1 ®— (2.4) 

where x a are local coordinates in M, ti,T2,T3 <E su(2) is a basis such that 

V{n,Tj) := -2Ti(TiTj) = 6ij, 

and r 1 , t 2 , t 3 is the dual basis. 

Einstein's theory of gravity is subject to constraints. In the standard ADM approach we have 
two constraints, namely the vector constraint generating the diffcomorphisms of M and the scalar 
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constraint generating dynamics, that is diffeomorphisms orthogonal to the Cauchy hypersurface 
M: 

C a {x) = Cf[x) + ■K{x) ( j>, a {x), (2.5) 

C(x) = C*(x) + ^^^ + ^\ x ) ( f >ta ( x ) ( j )!b ( x )^^ +V(4>)Vq(x), (2.6) 

where the terms Cf r and C gr involve the gravitational field variables q a b and p ab only. 
In LQG, the fields q a b and 

p ab in 

the constraints are expressed by the variables A\ and Pf , 
and we get an additional constraint - the Gauss constraint generating the "Yang-Mills" 1 gauge 
transformations of the fields (A, P): 

GHx) =d a P? + e ij k AiP%. (2.7) 

All the transformations generated by the vector, scalar and the Yang-Mills constraint are gauge 
transformations, because the constraints are of first class. 

In Section I and Section II the choice of the variables describing the gravitational part does 
not matter, so one can either use the ADM variables (q a b,P ah ) and the constraints (|2.5[ I2.6j) or, 
respectively, the Ashtekar-Barbero variables (A^ , P") and, the constraints (|2.5U2.6|l2.7|) . In Section 
III, the latter choice is necessary, because we will apply LQG. For the sake of the continuity, we 
will stick to the Ashtekar-Barbero variables, remembering that qab, p ab ,C gr and Cf r should be 
considered as functions of (A l a ,P°-). 

Each choice of the fields (A l al P°", <f>, it) defines a point in the phase space T. The solutions to 
the constraints form a constraint surface. We will also consider separately the phase space of 
gravitational degrees of freedom denoted by T gr , which by definition is set by the pairs (A l a , P"). 

By assuming that the vector and the scalar constraints are satisfied 

C(x) = 0, C a (x) = (2.8) 

we can solve the vector constraint in (|2.5j) for the gradient <fi M obtaining tfi ia = — and inserting 
this into the scalar constraint (|2.6[) . What we get, remembering (|2.8[) and solving the scalar 
constraint for it, is an expression for it 2 as a function of the geometry variables {A l a ,P°') and the 
potential V((f>) only, 

vr 2 = ^ (- (c GR + y/qv (<!>)) ± \j (c GR + V^m 2 - q ab c° R c^ . (2.9) 

The ambiguous sign ± in (|2.9[) defines different regions in the phase space T. In particular, only the 
choice of a plus sign includes the special case of a homogenous and isotropic geometry coupled to a 
scalar field. In the case of the minus sign specialized to cosmological spacetimcs, where each vector 
constraints vanishes identically, the expression for ir 2 above will just yield zero on the righthand 
side. 



B. A deparametrized model 

What we have done in the last section is solving the scalar constraint for the scalar field mo- 
mentum by using the vector constraint. Physically, this corresponds, as will be explained more in 
detail below, to choose the scalar field </> as our emergent time with respect to which the dynamics 
of the observables will be formulated. This calculation provides the relation between the standard 
real scalar field coupled to gravity, on the one hand, and the model we actually define below, on 
the other hand. 



Although we do not consider the Yang-Mills theory itself, the Ashtekar-Barbero variables are subject to the gauge 
transformations known from the Yang-Mills theory. 
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In our paper we will consider a model, that is denned by the vector constraint (|2.5|) . the Gauss 
constraint ()2.7j) and the following scalar constraint 

C'(x) = n(x)-h(x), (2.10) 

h := yj-y/qCsr + ^ (C^f - q^CFC? . (2.11) 

The scalar constraint C(x) has been rewritten using (|2.9[) . That theory is equivalent to the theory 
defined in the previous subsection in the case of no potential 

V((f>) = (2.12) 

and in the region of the phase space T such that '+' holds in (|2.9[) and 

tt>0. (2.13) 

Since the potential is set to zero in the model, <j> no longer occurs in the function h and the scalar 
constraints deparametrizes. Notice, that in the consequence of the constraints, in that region 

C gr < 0. (2.14) 

The deparametrized scalar constraints, being linear in the scalar field momentum, strongly Poisson 
commute 

{C'(x),C'(y)} = 0, (2.15) 

as a consequence of the following identity 

{h(x),h(y)} = (2.16) 

proved in [25| . A Dirac observable is the restriction to the constraint surface of a function / : T — s- R, 
such that 

{f,C a (x)} = {f,C'(x)} = {f,Gi(x)} = 0. (2.17) 

The vanishing of the first Poisson bracket means, that / is invariant with respect to the action 
of the local diffcomorphisms (that is all diffcomorphisms generated by the vector fields tangent to 
M) , the vanishing of the third Poisson bracket is equivalent to the Yang-Mills gauge invariance of 
/. The vanishing of the second Poisson bracket reads 

{f,n(x)} = {f,h(x)}. (2.18) 



III. QUANTUM CANONICAL GRAVITY COUPLED TO A SCALAR FIELD 

In this section we introduce a "formal" structure of our theory. Our goal, at this point, is to 
conclude what mathematical structures (Hilbert spaces, operators etc.) are needed to complete 
the quantization of the model. How to construct them using LQG will be explained in the next 
section. 

Assuming for the time being, that all Hilbert spaces and operators we need exist, and that they 
have the usual properties, we will now derive: 

• a general solution to the quantum constraints, 

• a general quantum Dirac observable, its classical interpretation and its physical evolution, 

• the Hilbert product between two solutions. 
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A. Quantum states and quantum fields 

The quantum states are complex valued functions 

OM) ^ *(<M), (3.1) 

where <j> and A are the scalar field and the Ashtckar-Barbcro connection defined on M in the 
previous section (henceforth, we will write A and P instead of A\ and P"). 
For a given representation the fields 4>, ir, A, P give rise to quantum operators 

fo)*(4>,A) = cf>{x)^{cj>,A), A) = IJtyfaA) (3.2) 

l o<p(x) 

A{{x)^A) = A{{x)*{4>,A) P*{x)*{4>,A) = l — S—UffaA) (3.3) 

These elementary operators are needed to define the operators corresponding to the classical con- 
straints and to define the quantum observables. 



B. The quantum constraints and their solutions 

We turn now to the quantum constraints and their solutions. The first step is defining the quan- 
tum counterparts of the classical constraints (|2.5|2.6I2.7|) . In LQG wc assume, that the quantum 
Gauss constraints corresponding to the classical expression in (|2.7|) still generate the "Yang-Mills" 
gauge transformations, hence their solutions are functions such that 

a~ x Aa + a' 1 da) = #(<M) (3-4) 

for every a : M -»SU(2). 

Similarly, we assume that the quantum vector constraints generate the local diffcomorphism 
transformations of the quantum states, and in the consequence, the quantum vector constraint 
carries over to the condition that ^> be invariant with respect to all local diffeomorphisms <p : M — > 
M, that is 

= 9(<j>,A). (3.5) 
The quantum deparametrized scalar constraint operator has the following form, 

C'{xyst = (tt{x) - h(x) \ (3.6) 

We use the equation (|2.11[) (which gives the expression for h as a functional of A and P) to quantize 
the second term in the parenthesis. Hcuristically we get 

h(x) = h(A,P)(x). 

Due to operator ordering aspects the definition of h is not unique and will be completed later in 
this paper. In order to avoid a quantum anomaly we must respect the classical symmetry in (|2.16[) 
also at the quantum level and must make sure, that 

[h(x),h(y)}=0, (3.7) 
(compare to (|2.16|0 . Given the quantum constraint operator (|3.6p . the constraint itself reads 

(jt(x) ~ h(xj\ * = 0. (3.8) 
To solve the quantum deparametrized scalar constraint, we write $ as 

* = ei/^WC*)^)^ (3.9) 
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with a new function ipi and insert it in (|3.8[) to obtain 
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*(</>, A) = ih(x)V((j),A). (3.10) 



54>(x) 

Due to the commutator in fl3.7[) . the constraint equation (|3.10p turns into 

ik)^ = °" (3 - 11} 

Hence, a general solution to (|3.10[) is 

#(<M) = e ^ d3x ^ x) ' h{x) ip{A). (3.12) 

Notice, that the exponentiated operator acting at t/j on the right hand side of ()3.9p is Yang-Mills 
gauge, and diffeomorphism invariant itself. Therefore: 

A general solution to the quantum vector, gauss and scalar constraints is every function \3.12\j . 
such that for every local diffeomorphism ip : M M, 

i>(<p*A) = HA), (3.13) 

and for every a : M -^-SU(2) 

iP(a- 1 Aa + a- 1 da) = ip(A). (3.14) 

In the remaining part of the article we will be using the abbreviation 

d 3 xfh := [ d 3 x<P(x)h(x). (3.15) 



C. Quantum Dirac observables 

A quantum Dirac observable is the restriction to the space of solutions to the quantum constraints 
of an operator O which satisfies the following two properties: 

• O is invariant under local diffeomorphism and Yang-Mills gauge transformations, 



[6,C'{x)\ = 0. (3.16) 

Following the ideas of the relational framework for observables 0-0] it is easy to construct a large 
family of Dirac observables. Let L be a linear operator which maps the functions A i-> ip(A) into 
functions A H> Lip(A). Consider an operator 

O(L) := ^fd 3 x4>hi e ~tfd 3 x4>h_ (3 -^ 

As required, the operator 0{L) commutes with the quantum version of the deparametrized scalar 
constraints, 

[0(L),C'(x)} = 0. (3.18) 

Moreover, the operator O(L) is Yang-Mills gauge and local diffeomorphism invariant provided the 
operator L is. 

Each of the operators O(L) defined by a Yang-Mills gauge, and diffeomorphism invariant operator 
L preserves the space of solutions to the constraints. Indeed, 

0(L)e^ d3xih i;(A) = e^ '^V' (3.19) 
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The operators (|3.19[) with the Yang-Mills gauge, and local diffcomorphism invariant operators L set 
a family (algebra, modulo the domains) of the Dirac observables. The total scalar field momentum 
J M d 3 xw(x) also defines one of the quantum Dirac observables (|3.19|) . namely 

0( [ d 3 xh{x)) = f d 3 xh(x). (3.20) 

JM JM 

The family of the Dirac observables p,19|) in fact contains all the quantum Dirac observables. 
To see that this is true, suppose an operator O satisfies the condition (|3.16[) at each x £ M . Let 
us write the operator in the following way 

Q = e ifd 3 x4,hj£ e ~tfd 3 x$h^ (3.21) 

where K is a priori arbitrary operator. The condition (|3.16j) with O substituted for the right hand 
side of (|3.21[) , takes the following form 

[K,n(x)} = 0. (3.22) 
The set of all the solutions K to (|3.22p is generated by the following solutions: (i) given any x <G M, 

K = 7r(x), 

and (ii) 

K = L, 

considered above, that is L which maps the functions A h-> tp(A) into functions A H> Li/j(A). 

The solutions of the type (ii) give rise exactly to the family of the quantum Dirac observables 
(|3.19[) we have introduced above. On the other hand, a solution of the type (i) gives rise to the 
following quantum Dirac observable 

However, we should keep in mind that what really defines a quantum Dirac observable is the 
restriction to the space of solutions to the quantum constraints. The restriction of (|3.23p is iden- 
tically zero. This shows that all the quantum Dirac observables are those defined by (|3.19p and 
diffeomorphism and Young-Mills invariant operator L. 



D. Classical interpretation of the Dirac observables 

Suppose, that a given operator L used to construct the Dirac observable 0{L) corresponds in 
the quantum theory to a classical function L defined on the gravitational phase space r gr , and that 
the support of L is contained in the set on which 



C* gr < 0. 



(3.24) 



To find a classical observable O(L) whose quantum counterpart is O(L), it is convenient to express 
the operator (|3.26j) in terms of a formal power series given by 



L, / d 3 x(j>h 



(3.25) 



(n) 



where [., .](„) denotes the iterated commutator defined by [L, J d 3 x(j>h]^ = L and [L, J d 3 xcj>h]^ 
[[L, J d 3 x<fih]( n _ 1 - j , J d 3 xcf>h]. The usual substitution [•, •] h->- •}, leads to a formal series 



oo 

°w = E / d3x ^ 



(») 



(3.26) 
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for a classical observable O(L). That series is very well known in the theory of relational observablcs 
0-H HH . To recall its meaning we first consider a slightly more general expression with <f> replaced 
by a point dependent parameter M 3 x t— > t(x), namely 

The * denotes the pullback, and the map 

a t : r gr -> r gr 

is defined by the hamiltonian flow [3 t : T — > T generated in the full phase space T by the constraints 
C'(x) with the parameters t(x). The action of the flow reads, 

Pt(A,P,<t>,T) = (a t (A,P),0-t,7r). 

Clearly 

MA,P,^tt) = ( a<t> (A,P),0,TT). (3.28) 
The value of O(L) at any point (A, P, <j>, it) is defined to be 

0(L)(A,P,^,7r) = L(a^(A,P)). (3.29) 

In conclusion, the quantum Dirac observable 0(L) corresponds to the classical function, that is 
also an observable, O(L), 

CHL) = O(L). (3.30) 

At this point a comment about the status of the operator O(L) is appropriate. It may happen, 
that given a point ( A, P, (p, it) in the classical phase space, the series Q3.26P is non-converging. In 
fact we encounter cases like that in the LQC models of homogeneous isotropic universe with positive 
cosmological constant |22|. However, still the operator O(L) is well defined as long as a self adjoint 
extension for the operator J cPxtfih is fixed, and therefore the unitary operator exp(i J d x (j)h) is 
well defined. Then, the quantum evolution just goes beyond the classical theory. That is exactly 
the reason, why we have chosen to define the Dirac observablcs directly in the quantum theory, 
and only interpret them in the classical theory as secondary objects. 



E. Dynamical evolution of the observables 

The Dirac observables we have defined are relational observables (often called "partial" 0,(1], 
section 1.2 of [J] ). For that class of observables one is able to define a non - vanishing evolution 
generated by a so called physical Hamiltonian, that will be introduced in the next section. The 
dynamics is defined with respect to an internal time given by the values, which that field (j> takes 
while being transformed along its gauge orbit. This can be seen in the construction of the quantity 
OiV) from a given function L by generalizing the choice of the evaluation point from (|3.28p to 

/3<t>-<t> (A,P,(l>,ir) = (a0_^ o (A,P)>o,7r), (3.31) 

where 4>o is an arbitrarily fixed function on M. We denote the resulting function defined on the 
phase space T by O^, (L), that is 

O^ (L)(A,P^tt) = L(a^ a (AP)). (3.32) 

For the function O,p (L) to be well defined, the flow f3 t : T — > T has to be well defined for 

t = 4> - 0o 

in the domain of the function L. 
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That classical construction leads to a corresponding quantum operator definition 

O 0O (L)*(0,A) = e iSd 3 x(4>{x)-M*))hx)i e -iSd 3 x{4>(x)-4> )Hx)^^^ A ^ (3 33-) 

where we used (j>^(4>, A) = 0^(0, A). This definition will not enlarge the class of the Dirac 
observables (|3.19[) . indeed 

0^(L) = O(L') (3.34) 

with 

Jj — £ -i f d 3 x4>a(x)h(x) j^ e i f d 3 xtj> h(x) (3.35) 

In this way, in the algebra of the (formal) solutions to the condition 

[0,C'(x)) = (3.36) 

we have defined an abelian group of automorphisms labelled by the functions 0o defined on M, 
namely 

O(L) h-> O^L). (3.37) 

If we want to restrict the automorphisms to the algebra of the quantum Dirac observables, we 
encounter an obstacle. Given a function <pQ, we want the operator (|3.35p to be diffeomorphism 
invariant for every diffeomorphism invariant operator L. For the operators h(x) that will be 
constructed from the LQG framework, that condition can be satisfied only for a constant function, 

</>o(x) = (f>o € K, for every x <E M . (3.38) 

The result is a 1-dimensional group of automorphisms of the algebra of the quantum Dirac observ- 
ables. The group encodes the dependence on the internal time of the algebra of the quantum Dirac 
observables. 

F. The physical Hamiltonian 



The dynamics is generated by the following equation 

d 

dcj) 

Where 



4>a(L) = -i[h phys ,O^ (L)} (3.39) 



/iphys := J d^xhix) (3.40) 

is usually called the physical hamiltonian for the reason that it is a non- vanishing Dirac observable 
generating true 'physical' evolution in contrast to the Hamiltonian constraint. 

The physical Hamiltonian will be an exact implementation of the heuristic formula 



i phys = / d 3 x\/ -VqCar + ^(Ctorja _ q^CTcT- (3.41) 

We remember however, that the operator will be applied to diffeomorphism invariant states (|3.5[) 
whereas the operator C^ r should generate the diffcomorphisms. Therefore, assuming the suitable 
choice of the ordering, the physical Hamiltonian acting on the diffeomorphism invariant functions 

V> is 



h phys iP(A) = J d 3 x^/-2^C9^(A), (3.42) 

where we also took into account (recall (12.141) ) , 

C gr < 0. (3.43) 

This result coincides with that of 1241. 
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G. The Hilbert product between the solutions: H p h ys 

Suppose we have a sesquilinear scalar product for the Yang-Mills gauge and local diffeomorphism 
invariant functions (or distributions) defined on the space of the Ashtekar-Barbero connections. 
Denote the product of the functions tp and tp' by 

(3.44) 

and the corresponding Hilbert space by Hdiff- 

We can use it to define the "physical" (that is respecting the dynamics) Hilbert product in the 
space of solutions p. 12ft : 

(e^Vle^'V) := (V#0- (3-45) 

V / phys 

The resulting Hilbert space "H p hy S is "physical" , and its elements are the physical states. 



H. Summary: the exact structures we need 



In summary, in order to construct the quantum model we will need: 

• the Hilbert space Hdiff of the Yang-Mills gauge and the local diffeomorphism invariant quan- 
tum states of geometry, 

• the operators in "Hdiff which admit a well understood geometric interpretation, 

• the physical Hamiltonian operator h p h ys defined in a suitable domain in Hdifi (which is not 
expected to be dense, because the heuristic formula for the operator involves the square roots 
of non definite expressions) . 

Given all that, the physical Hilbert space is unitarily isomorphic via 

with the domain of h p h ys in "Hdiff ■ 

Every observable O(L) (for simplicity let L be bounded) is the pullback by (|3.46[) of an operator 
L which preserves the completion of the domain of h p h ys - 

Finally, the emerged dynamical evolution (|3.39[) of the observables reads 

L{t) = e - ir Wi^W. (3.47) 
This is precisely the very well known Heisenberg picture evolution defined by the Hamiltonian 

^phys ■ 

Notice, that in fact, it is not necessary for L to preserve the domain of /i p h ys - Indeed, given any 
ip in that domain, the expectation value 

(iP\e- iT,l ^Le lT ' h ^^) = (e iT ' l ^ip \ Le"' h ^^) 

is well defined. This can be seen by using that it is equivalent to replace L by the operator 

V = PLP, (3.48) 

where P is the orthogonal projection onto the completion of the domain of ft. p hy S , and to considering 
the pullback of the Dirac observable O(L') together with its dynamics. 

This kind of structure will be necessary for the outcome. This is all we need to complete the 
quantization of a model of quantum gravity coupled to a scalar field. 

In the derivation of the operator corresponding to the physical hamiltonian /i p h vs , however, we 
will need yet more structure: 



11 



• the operator h p h ys should be defined by using the suitably defined operator valued distribution 

• the distribution should be self-adjoint, so that we can use the spectral decomposition to define 
the subspace 

y/q(jc)CP(x) < (3.49) 
and thereon the new operator valued distribution 



h{x) = y/-2y/q(x)C*{x), (3.50) 

• we should be able to verify the condition 

[h{x)Mv)\ = 0, (3.51) 

• and finally define 

Vy S = / d 3 xh(x). (3.52) 



Notice, that none of the operators ^/q(x)C ST (x) or h(x) can be defined within the Hilbert space 
Hdiff, because the x dependence manifestly breaks the diffeomorphism invariance. Therefore, the 
properties of the sclf-adjointncss require some extra Hilbert spaces, Hdif{,x, labelled by the points 
of M, whereas the commuting at different points can be defined only on a yet bigger Hilbert space. 

Remarkably, all the suitable structures can be constructed within the LQG framework, as we 
will explain in the next section. 



IV. APPLICATION OF LQG 
A. The Hilbert spaces 

1. The kinematical Hilbert space of quantum states of the geometry 

In LQG (we use the notation of 01), the kinematical Hilbert space of quantum states of the 
geometry is set by the so called cylindrical functions of the connection A. A cylindrical function is 
defined by a set a of finite curves e\, e„ in M and by a continues function / : SU(2)" — > C, in 
the following way 

1> a , f (A) = f(A( ei ),...,A(e n )) (4.1) 

where the symbol A(e) denotes the parallel transport along e defined by the connection A. The 
set Cyl of the cylindrical functions is a vector space, and an associative algebra. The space of the 
cylindrical functions Cyl is endowed with an integral 



ip a ,f i-> / rpa,f (4.2) 
used to define the sesquilinear scalar product 



(VwlVV,/')gr = / Tpaji>a'J>, (4.3) 

and defines (after the completion) the kinematical Hilbert space H. for the geometric operators in 
LQG. We assume in this paper that the manifold and the curves are piecewise analytic. Then, for 
every cylindrical function there exist curves a = {ei,...,e„} which form a graph embedded in M 
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(that is they are allowed to intersect only at the ends) such that the function is given by (14. ip . 
The curves ej are called edges of the given graph a. 2 
For a cylindrical function defined by a graph, we have 

VW = / d n gf( gi ,...,g n ), (4.4) 

JSU(2)™ 

where d n g is the Haar measure on SU(2) n . The geometric operators preserving the space Cyl are 

A(e)^ aJ (A) = A(e)^f(A( ei ),...,A(e n )) (4.5) 

and 

1 f 5 



Pi^aJ = 7T / TTjVw Vabcdx" A dx c . (4.6) 

There is an orthogonal decomposition of T~L into subspaces H.' a labelled by the embedded graphs 
a. To define it, denote first by (unprimcd) H. a C H the Hilbcrt subspacc spanned by the cylindrical 
functions ipa.f, with all the possible functions /. Those spaces, however, are too big to provide 
the orthogonal decomposition. Given a graph a, whenever a graph (3 can be obtained from the 
edges of a by glueing, or reversing the orientation or removing some of them, then Hp C H a . 
Therefore, define H' a C T-L a to be the orthogonal complement in H a of the subspacc spanned by 
those subspaces Hp. The decomposition is 

H = 0K, (4-7) 

a 

where a runs through the set of all the semianalytic embedded graphs in M. 



2. The Hilbert space of the diffeomorphism invariant states of the geometry 

Semianalytic diffeomorphisms Diff(M) of M preserve the space Cyl and act unitarily in the 
Hilbcrt space H just by the natural pullback of the Ashtekar-Barbero connections. Denote the 
action of ip eDiff(Af) by 

u^-.n^n. (4.8) 

To implement the construction of the quantum operator corresponding to the physical Hamiltonian, 
we will need two different Hilbert spaces: One of them includes states, that are invariant with 
respect to all (semi- analytic) local diffeomorphisms Diff(M) of M and the other one is the home 
of the states invariant with respect to the subgroup Diff(M, x), which preserves a given point 
x G M. (Later, we will also impose the Gauss constraint, that is the condition of Yang-Mills gauge 
invariance). Let Diff stands for either Diff(M) or Diff(M, x). The only Diff invariant direction in 
H is the constant function. However, since the group Diff is not compact, we expect the invariant 
states to be distributions on the space of the Ashtekar-Barbero connections, that is linear maps 

(*|:Cyl->C. (4.9) 

Whereas the space of all distributions seems to be too big, a suitable rigging map can be defined, 
which carries each ip £ Cyl into a Diff invariant distribution i/DiffCV 1 )- To recall the definition of this 
map, we need the orthogonal decomposition (|4.7j) . The map ?yDi£f is introduced for each subspace 
T-L' a individually. By the linearity, it extends to every cylindrical function. That is, the domain 
of the rigging map ?7Diff is Cylc H. The first step in the construction of the rigging map r?Diff , 



2 To be more precise, in what follows, an edge is either an oriented semianalytic imbedding of a circle in M, or 
a paramctrization free, oriented, finite curve defined by e : [0, 1] — > M such that either e is an imbedding, or 
e(0)=e(l). 
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is identification of the elements of W a that will be annihilated. Consider those diffcomorphisms 
if EDiff which map each edge of a into another edge modulo the orientation, and let us call them 
the symmetries of a and denote their group by Diff Q . The functions ip € H' a invariant with respect 
to Diff Q form a subspace denoted either by H! a inv in the Diff=Diff(M) case, or W a inv x , in the 
case of Diff=Diff(Af, x). The elements of 7i' a orthogonal to H' a inv are annihilated by the rigging 
map TjDiff- For ip € H' a inv , r?Diff(V ; ) is defined as follows 

TfrifftyO : i>" ^ ( U ^W)- (4.10) 

[4>\ GDiff /Diff Q 

Note, that if ip" G , then the right hand side is zero if a" is not Diff equivalent to a, and in 
the case there is <fi" GDiff such that 

<j>"{a) = a", (4.11) 
the only possibly non-zero term in the sum in (|4.10j) is 

VdsbW : V>" (EVW)- (4-12) 

Since every cylindrical function is a finite sum of elements of the Hilbcrt spaces H' a , ?7Diff (ip) is 
defined in Cyl. For the same reason, the map 

^i-HftjjfftyO (4.13) 

extends by the finite linearity to Cyl. 

With the rigging map ryDiff we define not only the vector space of the Diff invariant states to be 
the image 77Difr(Cyl), but also the sesquilinear product 

(?7Difr(^)|f7DifT(^'))Diff := (j7Diff(V0> V0 (4-14) 
In this way we have defined a Hilbert space "HDiff ■ The map ?7Diff defines a natural isometry 

n DiS ee 0K, Diff (4.i5) 

M 

where [a] runs through the set of the Diff classes of the graphs embedded in M. Recall that Diff 
= Diff (M), Diff(M, x). Therefore, we have defined two types of the Hilbert spaces: the Hilbcrt 
^Diff(M) an d, respectively, per each point x £ M. the Hilbcrt space %-ois(M,x)- 

3. The Hilbert spaces of the Yang-Mills gauge and diffeomorphism invariant states of the geometry 

Imposing the Gauss constraint is yet easier, than requiring diffeomorphism invariance, and could 
be equivalently done, either before, or after solving the diffeomorphism constraint. The group of 
unitary transformations of % given by the Yang-Mills gauge transformations is compact. Hence all 
solutions to the Gauss constraint in % arc invariant elements of H. (as opposed to non-normalizable 
states, distributions). Moreover, the group of the Yang-Mills gauge transformations (|3.4j) preserves 
each of the subspaces "H' a . For every Yang-Mills gauge invariant ip sCyl, the Diff invariant distri- 
bution 

musty) G Hdib (4.16) 
is also insensitive to gauge transformations of ip" GCyl. Namely, the number 

mxstyW) 

is invariant. The converse is also true: If ??Diff (ip)(tp") is invariant with respect to the Yang-Mills 
gauge transformations of ip" ', than ip is Yang-Mills gauge invariant. 
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In conclusion, the Yang-Mills gauge and diffcomorphism invariant distributions on the space of 
the Ashtekar-Barbero connections we want to use to construct the Hilbert space HdiS of section 
III.G, are the distributions 

%>iff(M)(V0 (4-17) 

obtained from the Yang-Mills gauge invariant cylindrical functions ip. Denote their Hilbert space 
by Hdiff- By construction 

Hdis C Huiff(M)- (4-18) 

For the introduction of the physical Hamiltonian we will also use the Hilbert space Wdm.x 
obtained by replacing in the construction of the Hilbert space "Hdiff the group Diff(M) by the 
group Diff(M,x). 

B. The operators 

1. The Dirac observables 



From the previous subsection we already have the LQG candidate for the Hilbert space Hdis 
of the Yang-Mills gauge invariant and diffcomorphism invariant quantum states of geometry. As 
we already know from Section III.G, from a suitable subspace of this space we will construct 
the "physical" Hilbert space of solutions to all the constraints of the model we are considering. 
Secondly, in the Hilbert space "Hdiff we will need the operators representing the geometry of the 
initial data defined on M , from which we will construct the Dirac observables. 

Let us begin with this second task, because it is easier. We assume below, that the operators 
we consider in the Hilbert space Ti, as the domain have the vector subspace Cyl of the cylindrical 
functions. Every Yang-Mills gauge and Diff(M) symmetric operator L defined in the kinematical 
Hilbert space "H, defines naturally by the duality a symmetric operator L in Udm, 

(^toiff(M)W>)>"> := (?7Diff(M)(V0,^"} = (rtom(M)(Lip),tp") (4.19) 

where the bracket is the action of a distribution (a first entry) into a given cylindrical function xp", 
that is, we could phrase it in a simpler way 

Lr)r>is(M)(^) = Vnitt{M)(Lil>). (4.20) 

An excellent example of a Yang-Mills gauge and diffcomorphism invariant operator in Ji available 
in the literature [H^HH is the volume of the underlying manifold M operator 

Vm = J d x jq(x). (4.21) 

Another example we manage to construct might be any quantum operator representing the integral 
of a scalar constructed from the intrinsic or extrinsic curvature. 

In the kinematical Hilbert space H, there is also a well defined operator valued distribution 

jq(x) = Yl S(x,x?)jq xn (4.22) 
x'eai 

where each of the operators y/q , is Diff(M, x r ) invariant. The uncountable sum on the right hand 
side is well defined, because for every smearing function F : M — > R, and a cylindrical function 
ipaj, we have 

r n 
d 3 xF(x)^q(x)ip aJ = J2 F ( v M Vl i><*,f, (4.23) 
^ i=i 

where V\, ...,v n are the vertices of the graph a. Via (|4.20[) . for every x' £ M, the operator y/q , 
defines an operator yfq x , in Jidiff^'- Morally, y/q(x) is also Diff(M, x) invariant for every given 
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x g M, therefore (|4.20[) should also be somehow generalized to this case. Indeed, (see [Hj]) the 
suitable generalization is natural and provides in this case a distribution 

Jq{x) = S(x,x')^ x „ (4.24) 

x'eM 

which makes sense due to the fact that all the Hilbert spaces ridis,x are embedded in the single 
vector space Cyl*. 

There is one more technical remark in order at this point. Consider two operator valued distri- 
butions in 'H, of the form 

A(x) = J2 S ( x , x ')Ax>, B(x) = J2 S ( x , x ')Bx> (4.25) 
x'eM x'eM 

each of which satisfies the property (|4.23|) . A natural regularization by smearing leads to a new 
operator valued distribution 



A(x)B(x) = J2 S(x,x')\/S(A X ,B X ,) (4.26) 
x'eM 

which also has the property (|4.23|) . where S stands for a symmetric product of the operators, and 
the domain of the resulting operator is restricted by the positivity of S{A X B X ) requirement. The 
regularization consist in the smearing 

Mx) = J d 3 A(y)6 £ (y,x), B e (x) = J d 3 B(y)S e (y,x) (4.27) 

with a smearing function whose support goes uniformly to x = y as e — > 0, and which goes to the 
Dirac 6(x,y). The key trick is an observation that for every fixed graph a, for sufficiently small e 

n 

A e (x)B e (x)ij a j = ^(^(.t,^)) 2 ^^^,/, (4.28) 

7=1 

for any cylindrical function ipaj-, and more over, the sum on the right hand side contains at most 
one non zero element. Due to the latter property 



A t (x)B € (x) = J2 6 *( x > v i)\/ A VI BvA a ,f (4.29) 



i=i 

provided the square root is well defined itself. Finally, 



d 3 xF{x)yjA e {x)B e {x)ij a j ->■ ^F(vi)^A VI B Vl ip aJ . (4.30) 



i=i 



2. The quantum scalar constraint and the physical Hamiltonian 

As we remember, our first task we can finally turn to now, is a construction of the physical 
Hamiltonian operator 

W = J d 3 x^-2^(xjc&(x) (4.31) 

defined in Hdis- 

A quantum scalar constraint C gr was defined in [30l ]. its properties and possible generalizations 
were studied in 0, HH . We will be using here the formulation of the scalar constraint of [|| . In order 
to use it for our current construction, we will need a new element. Thus far, the scalar constraint was 
used either as smeared against arbitrary laps function J d 3 xN(x)C(x), or, as the master constraint 

J d 3 Xyfq{x)~ 1 C{x)C^ {x), or as a physical Hamiltonian defined after deparametrization with respect 
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to 4 scalar fields. The smeared scalar constraint maps a domain in "Hdiff into Cyl*, there is no 
sense in which it could be symmetric or sclf-adjoint. The master constraint, on the other hand, as 
well as the physical Hamiltonian after the 4-fold deparamctrization, respectively, is defined in the 
kinematical Hilbert space W as a graph preserving operator. The current case, is a new one, we 



will need an operator y — 2 s /q(x)C sr (x) defined in Haiti- 

The quantum scalar constraint presented in Q takes the following form, 

/ d 3 xN(x)C(x) = C x , (4.32) 
J xeM 

where each of the operators C x maps its domain contained in 'Hdiff into 1-Ldi«,x- However, as it 
follows from [32jj, it naturally defines an operator in the corresponding Hilbert space "Hdiff, The 
advantage is, that only now we can require the symmetry (self-adjointness) of those operators. As 
defined in the operators C x come out non-symmetric. The minor improvement, but necessary 
for our current work, consists in replacing them by symmetric operators 

Cf = \(C X +Cl). (4.33) 

and choosing an essentially self-adjoint extension that may be non-unique. Then, the quantum 
scalar constraint operator we will use for the physical Hamiltonian takes the following form 

On the other hand we have already considered above the volume density quantum operator written 
in the similar form, 

Vffj = S(x,x')^ x ,. (4.35) 

x'eM 

At this point, we are in the position to define the operator 



-2y/q(x)C&(x) (4.36) 

A regularization in H similar to the one discussed above, gives (modulo the symmetrization of the 
product of the operators yfq^i and C^) 



-2^q(x)C&(x) = Y S(x,x')\l-2^ 1/2 C^^ 1/2 =: h(x). (4.37) 

x'eM 

However, the operator is well defined only in the subspacc of "Hdiff. x corresponding to the positive 

-——1/2 . _ 1/2 

part of the spectrum of ^JqZ C| r \/fe • To formulate that condition we need to choose a self- 
adjoint extension of the operator in the case it is not unique. Denote the resulting subspacc of 
"Hdiff, Z by %dxS,x+- There is a natural averaging map 

VM ■ HdiB,x -> UdiS (4.38) 

VDiS(M,x)W l-> ?7Diff(A/)(^)- (4.39) 

The domain of the physical hamiltonian is 

"Hphys = TlM{%iiB,x+), (4.40) 
and the formula for physical Hamiltonian reads 



^phys 



f , - v / —1/2. -1/2 

/ d 3 xh(x) = V^Vfe Cfjqx ■ (4.41) 

J XEM 
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We remember the anomaly free condition f|3.51[) that should be satisfied by our construction. In 
(32l | an extension of the Hilbert space Hphys is introduced in which the smeared scalar constraint 
operators 

C sr (N) = [ d s xN(x)C SI {x) (4.42) 

J M 

are defined and their products 

C sr (N)C ST {N') (4.43) 

are contained. It follows from the results of [1H that the smeared constraint operators commute 

[C SI {N),C sr {N')} = (4.44) 

on a large subspace of the enlarged vector space. It justifies our conjecture, that the condition 

[h{x)Mv)\ = (4.45) 

is a restriction on the ambiguities in the definition of the operators h(x), that is on the loop 
assignment 0, H| and the self-adjoint extensions. 



V. CONCLUDING REMARKS, OUTLOOK 



We have another quantum model model of gravity involving all the degrees of freedom. The 
model discussed here assumes a vanishing potential for the scalar field that becomes the internal 
time for the Dirac observables. Neglecting this requirement has the effect that the physical 
Hamiltonian depends on the internal time </> as can be seen in equation (|2.9[) . Non - conservative 
Hamiltonians usually increase the intricacy as far as the technical perspective is concerned. 
Likewise if we use for instance Standard Model matter instead of a scalar field the system will also 
not deparametrize anymore. Hence, all the technical simplifications due to deparamctrization used 
in this work are not available any longer. A discussion about which kind of matter Lagrangians 
induce a deparametrization for General Relativity can be found in [l2j . 

The quantization of this model is complete and every necessary element exists within the 
framework of LQG. However, there are still ambiguities though, present in the LQG definition 
of the quantum scalar constraint operator due to its non-polynomial structure. The only way 
to understand them and their possible physical meaning is to start applying the model. Before 
explaining what the model discussed in this work is good for, let us compare it briefly to the first 
model that was completed by Giesel and Thiemann. 



A. Comparison with the Brown-Kuchar model applied to LQG 

The Brown-Kuchar (BK-) model 0] considers four scalar fields that have the properties of dust 
and become a dynamically coupled observer, with respect to which the dynamics of the remaining 
degrees of freedom is formulated. This model was used by Giesel and Thiemann [23| and a reduced 
phase space of gravity coupled to dust was derived. For this purpose the BK-model needed to be 
extended since the reduced phase space requires also the construction of (classical) Dirac observ- 
ables with respect to the scalar constraint. The original BK-model is rather the counter part of 
what is done in this paper because there the vector constraint was reduced classically, whereas for 
the scalar constraint a quantum condition was formulated. 

In the reduced phase space quantization procedure discussed in [23| both, the scalar as well as 
the diffeomorphism constraint, are reduced classically. The Gauss constraint is, as in this paper, 
solved at the quantum level. This yields to an algebra of observables describing the classical phys- 
ical phase space. Due to the deparametrization of the scalar constraints, this algebra turns out to 
be isomorphic to the kinematical one. In contrast to what is done in this paper, a quantization of 
the observable algebra accesses directly the physical Hilbert space (once also the Gauss constraint 
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is satisfied). Since the kinematical algebra is isomorphic to the physical one, in [23[ the standard 
kinematical representation of LQG can also be used for the physical Hilbert space 'Hphys- Similar to 
the work in this paper, the generator of the physical dynamics, the so called physical Hamiltonian 
ft-phys, is invariant under local diffeomorphisms. In the reduced approach this leads to the require- 
ment, that in order to avoid a quantum anomaly, the operator needs to be invariant under local 
diffeomorphisms too. As shown in (34j operators being invariant under local diffeomorphisms and 
defined in the standard (kinematical) LQG representation cannot be graph-changing. This means, 
that they need to preserve the graph they are acting on, yielding the condition, that the LQG con- 
straint operators 0, entering the physical Hamiltonian h p hy S m [HI need to be quantized in a 
graph-preserving way. As we explained above, LQG is glued from the Hilbert spaces corresponding 
to all possible graphs. The original LQG scalar constraint operator does not preserve those graph 
Hilbert spaces. In the model of [23| the physical Hamiltonian must preserve each graph Hilbert 
space. In the consequence, the constraint operator has to be suitably redefined in [23( when the 
standard (kinematical) LQG representation is used for 'Hphys- The paper [23| also discusses the 
quantization of the reduced model in the framework of Algebraic Quantum Gravity [HI], where 
a different representation is used, namely von Neumann's infinite tensor product representation. 
The quantum dynamics is not defined on embedded graphs but on abstract ones, carrying only 
combinatorial information. In this framework only the infinite combinatorial graph, that the the- 
ory is defined on and that acts like an abstract lattice, needs to be preserved by /iphy S) whereas 
any possible subgraph of this does not. In the case of the model presented in this paper here, the 
graph Hilbert spaces are not preserved and they evolve in the emergent time. 

B. Application of this model 

Our model can be used to verify the properties of quantum space-time we expect after learning 
the lessons from LQC and QFT in curved spacetime. 

In the LQC models of the homogeneous masslcss scalar field coupled to gravity, Big Bang turns 
out to be replaced by Big Bounce, as the result of the quantum gravity effects. Now, with our model, 
we can consider the same system of fields from the point of view of the full theory, without the 
symmetry reduction. Similarly, we can also consider the quantum gravitational collapse, quantum 
black holes, theory entropy. All those cases are manageable within our model, and the only difficulty 
is of technical nature. Also the Hawking radiation and black hole evaporation process expected 
from the theory of quantum fields on the classical black hole background are in the range of our 
model. The next step to obtain progress in this direction is the construction of scmiclassical states 
for full LQG, which are preserved under quantum dynamics generated by the physical Hamiltonian 
on appropriate time scales. 

In conclusion, our paper opens the door to understanding the properties of quantum spacetime 
from the point of view of the full quantum gravity. 
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